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The electronic and magnetic properties of monoclinic double perovskite Sr2CeIr06 were exam¬ 
ined based on both experiments and first-principles density functional theory calculations. From the 
calculations we conclude that low-spin-state Ir 4+ (5d 5 , S=-|) shows t 2 ^ band derived anti-ferro type 
orbital ordering implying alternating occupations of d yz and d xz orbitals at the two symmetrically 
independent Ir sites. The experimentally determined Jahn-Teller type distorted monoclinic struc¬ 
ture is consistent with the proposed orbital ordering picture. Surprisingly, the Ir-5d orbital magnetic 
moment was found to be ~ 1.3 times larger than the spin magnetic moment. The experimentally ob¬ 
served AFM-insulating states are consistent with the calculations. Both electron-electron correlation 
and spin-orbit coupling (SOC) are required to drive the experimentally observed AFM-insulating 
ground state. This single active site double perovskite provides a rare platform with a prototype 
geometrically frustrated fee lattice where among the different degrees of freedom (i.e spin, orbital, 
and lattice), spin-orbit interaction and Coulomb correlation energy scales compete and interact with 
each other. 

PACS numbers: 71.15.Mb, 71.20.Be, 75.25.Dk, 71.70.Ej 


Transition metal oxides (TMOs) with orbitally degen¬ 
erate d electrons exhibit rich electronic and magnetic 
properties arising from the interplay of spin, orbital, 
charge, and lattice degrees of freedom. When the num¬ 
ber of d electrons (or holes) in a transition metal is 
smaller than the degeneracy of d orbitals, the d elec¬ 
trons have the freedom to occupy any of them. In 
such cases, the occupation of d orbitals at one site is 
often correlated with the occupation at the next site. 
Therefore, analogous to spin ordering, the orbital de¬ 
grees of freedom can be locked and particular orbitals 
are occupied at each site, which is called orbital order¬ 
ing. This orbital ordering, generally results from the in¬ 
terplay between structural distortion and Coulomb cor¬ 
relation of electrons. Orbital degrees of freedom often 
couple with spin degrees of freedom [1] and give rise to 
anisotropic electronic and magnetic interactions. Due 
to strong correlations between different degrees of free¬ 
dom, orbital ordering spontaneously drive various ex¬ 
otic phenomena such as reduced dimensionality]!!, struc¬ 
tural phase transition [2], exotic magnetism Q, Peierls 
transitions Q, charge ordering [H| and orbital ordering in¬ 
duced ferroelectricity[(jj. 

The possibility of orbital ordering has been dis¬ 
cussed for doubly degenerate e g orbitals of perovskite 
manganitesQ, cuprates El [j| and for the degenerate t 2 ^ 
orbitals of vanadates (§1, [lo|. Most studies investigating 
this aspect have focused on 3d TMOs; much less is known 
about whether related 4 d and 5 d oxides can exhibits sim¬ 
ilar behaviour. Recently, 4 d and 5 d TMOs have gener¬ 
ated significant excitement because of the observation of 
unexpected orbital ordering [111 and localized transport 
properties and magnetism (l2l. Il3j in ruthenate and iri- 
date compounds, respectively. Due to the large spatial 


extension of 5 d orbitals, the effects of electron correlation 
are weaker in 5 d oxides than in 3d oxides, while the trend 
is opposite for the effects of spin-orbit coupling (SOC). 
Therefore, 5d oxides should exhibit a different balance 
among the effect of spin, orbital, and charge degrees of 
freedom. Although the importance of SOC in magnetic 
materials has been well known, the interplay between 
SOC and electronic Coulomb interactions in magnetic 
materials remains a largely unexplored field. The rea¬ 
son for this seems to be the mutually exclusive domain 
involved in their acting. It is only recently that this sep¬ 
aration has been auestioned[12Hl5]. 

The 5d double perovskite family is one of the most 
widely investigated oxides for studies exploring exotic 
phenomena such as magneto-o ptic effect[l6], high-T c 
ferrimagn ets 171. half metallicity[18[, peculiar AFM spin 
structures [191421| . as well as strong SOC driven phe¬ 
nomena such as spin liquid [0, [23] , topological Mott 
insulator (lii and the Weyl semimetalpjj. In the con¬ 
text of orbital ordering, 5d based systems have been less 
investigated, because generally for 5d elements such as 
Os, Ir, the SOC energy scale considerably exceeds the 
structural distortion energy scale. Strong SOC destroys 
orbitally degenerate electronic states required for orbital 
ordering. In this letter we report a particularly striking 
example of the interplay among SOC, Coulomb interac¬ 
tions and crystal field splitting in the double-perovskite 
Sr 2 CeIr06, which shows a rare type of orbital order¬ 
ing originating from the Ir-t 2 ^ states. Using transport 
and magnetization measurements in addition to density- 
functional theory based first-principles calculations we 
explored this interesting aspect in the Sr 2 CeIr06 dou¬ 
ble perovskite compound. We propose a rare case of t 2 ^ 
orbitals derived anti-ferro type orbital ordering in this 
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FIG. 1: (Color online)(a) Crystal structure of Sr2CeIr06- 
Green, yellow, blue and red spheres represent Sr, Ce, Ir, and 
O atoms respectively, (b) and (c) Bond lengths marked oc¬ 
tahedral entities of Ir06 and Ce06 respectively, (d) Three- 
dimensional magnetization density plot showing the orbital 
ordering. The orbital characters d yz (pink surface) and d xz 
(green surface) at each Ir site are marked. The isovalue was 
chosen to be 70 e~/A 3 . The dashed lines designate the orbital 
chains in the structure. 


material at the Ir site with the half filled (5 d 5 ) config¬ 
uration while the Ce remains in the inert [Xe] configu¬ 
ration. Our proposed t 2 g band derived anti-ferro orbital 
ordering constitute a strong contrast to the most widely 
known anti-ferro orbital ordering in KCuF 3 [l] originated 
from the two e g orbitals. Moreover only a few frustrated 
magnetic materials have been investigated in the context 
of orbital ordering, such as NaNiOoj^. LiNi02 27, 2^] 
and spinel ZnV204[29[. In this context S^CelrCU double 


perovskite is a rare example of a prototype geometrically 
frustrated fee lattice where orbital ordering plays a cru¬ 
cial role. Our calculations revealed that the experimen¬ 
tally observed AFM insulating ground state is driven by 
strong SOC in addition to small onsite electron-electron 
Coulomb correlation at the Ir site. These particular be¬ 
haviours make this system a potential candidate for in¬ 
vestigating the interplay among the multiple degrees of 
freedom, such as spin, orbital, lattice, and multiple en¬ 
ergy scales e.g. SOC and Coulomb correlations, that in¬ 
teract and compete with each others. 

Sr 2 CeIr 06 crystallized in a monoclinic crystal struc¬ 
ture with space group P2i/n, with the lattice param¬ 
eters a =5.8255 A, b = 5.8445 A, c = 8.2435 A and 
/3= 90.17C these values are very similar to the reported 
structure [30|. Due to lack of sensitiveness of positions of 
light atoms such as oxygen in X-ray diffractions based 
structural analysis, the theoretically optimized structure 
of Sr 2 CeIr 06 was obtained by relaxing the atomic po¬ 


sitions of all atoms while keeping the lattice parameters 
fixed at the experimentally determined values. The struc¬ 
ture consists of alternating corner sharing Ir06 and Ce06 
octahedra, in all directions, as shown in Fig. 1(a). The 
Sr atoms are situated at the void positions between the 
two types of octahedra. Both Ir06 and Ce06 are dis¬ 
torted, the respective six metal-oxygens bond lengths are 
grouped into three values as shown in Figs. 1 (b) and (c) 
respectively. Interestingly the in-plane Ir-0 bond lengths 
of two neighbouring IrC>6 octahedra are of opposite type, 
i.e the short in-plane Ir-0 bond of one octahedron cor¬ 
responds to the long in-plane bond of the neighbouring 
octahedron. Though the positions of the oxygen atoms 
shifts in the optimized structure compared to the experi¬ 
mental values, the average Ir-0 and Ce-0 bond lengths of 
the theoretically optimized structure (2.011 A and 2.237 
A respectively) are in good agreement with the experi¬ 
mental average Ir-0 (2.009 A) and Ce-0 (2.202 A) bond 
lengths. 

The magnetic susceptibilities of S^CelrCU were mea¬ 
sured using a SQUID magnetometer (Quantum Design, 
MPMS MultiVu Application) in the temperature range 
from 2 to 300 K in varying field’s up to 7 T, see Fig. 2. 
Inverse magnetic susceptibility data indicated AFM or¬ 
dering at about 21 K. Fitting the high-temperature sus¬ 
ceptibility to the Curie-Weiss law [x = j&q], gave fi e ff = 
1.59 ii b with a Weiss constant of 0= -108 K, which are 
very similar to the previously recorded values |30[. A neg¬ 
ative value of 0 implies dominant AFM interactions in 
the system. The observed effective magnetic moment in¬ 
dicates, that Ir 4+ (d 5 ) is in S=| low-spin state and that 
is somewhat smaller than the expected spin-only 
contribution [^/S(S + 1)] for the S=| state. The inset 
of Fig. 2 shows the exponential decay of the resistiv¬ 
ity as a function of temperature measured in a 0.01 T 
field with the Quantum Design PPMS MultiVu Applica¬ 
tion. Fitting the high temperature resistivity data with 

E a 

the semiconducting hopping mechanism, p = poe 2K B T , 
resulted in a bulk band gap (E g ) of ~ 0.3 eV, indicating 
small gap insulating behaviour of this material. 

The calculated ferromagnetic (FM), GGA+U density 
of states (DOS) projected on to Ir-5d states is shown in 
the top left panel of Fig. 3. The empty Sr states lay far 
away from the Fermi level (E/), which is consistent with 
the nominal Sr 2+ valence state. At the E f only Ir-5d 
states contribute whereas the Ce-4/states are empty and 
situated far above the E/, although Ir-5d states show 
significant mixing with the Ce-4/ and 0-2 p states. We 
found that in the absence of on-site Coulomb correlation 
U, the band structure shows a metallic character. The 
incorporation of small Coulomb correlation U e ff > 2 eV 
at the Ir site, drove the insulating band structure. To 
clarify the effect of Coulomb correlation U on the elec¬ 
tronic structure we systematically varied U\f f from 0 to 
5 eV, as shown in Fig. S2 (in the supplemental mate- 
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FIG. 2: (Color online)Temperature dependence of the in¬ 
verse magnetic susceptibility of Sr 2 CeIr 06 measured at field 
strengths of 3.5 T and 7 T. The inset shows the temperature 
dependence of the resistivity measured in a 0.01 T field. The 
green (light) dots represent the experimental data while the 
solid line represents the fitting with exponential variation. 


rial). The gap, and magnetic moments [0.63 (Ib ( 

= 0 eV) to 0.86 /jLb (f/g// = 5 eV)] at Ir sites, increased 
with increasing £/, while Ce stayed in the nonmagnetic 
state with zero magnetic moment. The orbital occupan¬ 
cies and spin states remained unchanged across the entire 
applied U range. Comparing our calculated band struc¬ 
ture to the resistivity derived band gap, we choose to 
use U T e r ff = 3 eV (which, also is the commonly used U 
value for Ir in the literature |32j), for the following calcu¬ 
lations. Because of the octahedral environment the five 
d states were broadly split into t 2 g (d xy , d yz , d xz ) and e g 
(d 3z 2 _ r 2 , d x 2_ y 2 ) states as marked in the top left panel 
of Fig. 3. The Ir-t 2 P states were completely filled in the 
majority spin channel and partially filled in the minority 
spin channel, whereas the Ir-e^ states were completely 
empty in both spin channels. These findings in addition 
to the calculated magnetic moment suggest that Ce is 
in the 4+ (Ce 4+ :[Xe]) valence electronic state with zero 
magnetic moment in the inert configuration, and that Ir 
is in the 4+ (5 d 5 ) low spin valence electronic state with 
S =|. These findings are consistent with the experimen¬ 
tal magnetization measurements and previous electronic 
structure results [Hj. However in the previous studies 
one failed to drive the insulating band structure without 
inclusion SOC. 

Because of the octahedral distortion, the t 2 g and e g 
states are no longer degenerate. Contraction of Ir06 oc- 
tahedra along the c axis splits the triply degenerate t 2 g 
states into a lower d xy band (ai g symmetry), and upper 
d yz ~ d xz bands (e£ symmetry), as shown in the level 
diagram plot in the right top of Fig. 3. The d xy band is 
filled in both spin channels and this is electronically in¬ 
active. The two e^ bands are filled by two majority spins 
and one minority spin, which results S=-| state. However 


the single minority spin has the freedom to occupy any 
orbital between the double degenerate d yz — d xz bands. 
Therefore this minority spin experiences orbital degener¬ 
acy and drive a definite orbital ordering pattern in the 
lattice. To clarify this orbital degeneracy we plotted the 
orbital projected DOS for two Ir sites, Ir-1 (0.5, 0.0, 0.5) 
and Ir-2 (0.0, 0.5, 0.0), in the unit cell, as shown in the 
left middle and lower panels of Fig. 3 respectively. The 
plots show that the d yz state of the Ir-1 site is equiva¬ 
lent to the d xz state of the Ir-2 site, i.e. at the Ir-1 site 
the d xz state is occupied while for the Ir-2 site, the d yz 
state is populated by the single minority spin in the ej 
bands, which is clearly shown in the energy level diagram 
of the respective right panels. To visualize the pattern 
of orbital ordering we plotted the magnetization density 
within GGA+ U( U T e r ff = 3 eV) as shown in Fig. 1(d). This 
figure shows that crystallographically two non-equivalent 
Ir sites, i.e Ir-1 and Ir-2 in the unit cell are decorated with 
the alternating d yz and d xz type orbitals and that these 
alternating orbital chains run along the body diagonal to 
the monoclinic unit cell. Because in this case the Ir sites 
prefer two different types of orbitals, we described this 
type of orbital ordering as anti-ferro type orbital order¬ 
ing. Interesting point to be noted is that, aforementioned 
anti-ferro orbital ordering consisting of two t 2 g orbitals 
is in contrast to the most widely known classic exam¬ 
ple of anti-ferro orbital ordering as found in KCuFsJll, 
where alternating sites of Cu are occupied by the two e g 
orbitals. Therefore Sr2CeIr06 is a very unique case. 

The resulted orbital ordering could be understood as 
an effect of the cooperative Jahn-Teller distortion associ¬ 
ated to the Ir06 octahedra. The individual IrC>6 octahe- 
dra are distorted as reflected by the three different types 
of Ir-0 bond lengths in each octaherdon. If we concen¬ 
trate only on the ab plane Ir-0 bond lengths, theses are 
grouped into two Ir-0 bond lengths as short 1.986 A and 
long 2.033 A. However the crystallographic orientation 
of short and long bond lengths for the two Ir sites (i.e Ir- 
1 and Ir-2) are opposite. For Ir-1 the short bond (1.986 
A) is along the crystallographic b axis and the long bond 
(2.033 A) along the crystallographic a axis, while for Ir-2, 
short bond is along the crystallographic a axis and long 
bond along the crystallographic b axis, i.e short bond of 
the Ir-1 octahedron corresponds to the long bond of the 
Ir-2 octahedron, see right middle and bottom panel of 
Fig. 3. Therefore this cooperative Jahn-Teller distortion 
in the Ir06 octahedra, in terms of crystallographic ori¬ 
entation of Ir-0 bond lengths, induces the occupation of 
two different orbitals in the two neighbouring distorted 
Ir octahedra. 

To understand the role of SOC in this orbit ally ac¬ 
tive material, we performed GGA+H+SOC (^4// = 1~5 
eV) calculations. Surprisingly we found that GGA+SOC 
without onsite [/, could not drive the insulating nature 
of the band structure. We indeed required a small on¬ 
site correlation (— 1 e V) along with SOC to open 
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FIG. 3: (Color online) GGA +U(U{ r ff = 3 eV) FM spin- 
polarized Ir-5d orbital projected DOS. The two channels of 
the all panels represent majority and minority spin channels. 
The top panel represents the total Ir -d DOS in the unit cell 
while middle and lower panels represent the orbital projected 
DOS separately for the two Ir sites in the unit cell. The Fermi 
level is marked at zero on the energy scale. The right panels 
represent the energy level diagram corresponding to the left 
panel DOS. The individual Ir06 octahera in the right panels 
are shown to clarify the bond length orientation in the crys¬ 
tallographic axes. In plane short (1.986 A) and long (2.033 
A) bonds were marked by green and orange respectively. The 
arrow indicates the contraction and elongation of bond length 
along crystallographic axes. 


an insulating gap in the band structure. Because of the 
orbital degeneracy of Ir, we found that orbital magnetic 
moment to be larger than the spin magnetic moment at 
the Ir site, although the value of the orbital and spin mag¬ 
netic moments depended sensitively on chosen U values. 
For example, as we increased from a small value of U T e r ff 
= 1 eV \/j J spin (ii orb ' ltal )=0.25 (0.32) /xg/Ir] to moderate 
value of Ui r ff = 3 eV ^orbitai j = 0.33 ( 0 . 43 ) 

to a large value of ifjff = 5 eV [/i spm (/i or6ztaZ )=0.38 
(0.53) iib/ Ir], the spin, orbital moment and the ratio 
of ^orbital/^spin [ L2 g ( for ^ /= 1 e V), 1.34 (for U{ r ff = 

3 eV), 1.39 (for =5 eV)] increased. We note that 
the orbital magnetic moments point in the same direc¬ 
tion as the spin magnetic moment which is expected for 
the more than half filled t 2 p occupancies at the Ir sites. 
Similarly large orbital magnetic moment and a similar 
fjorbital j^svin ra ti 0 are reported in the literature in the 
context of an other Ir based perovskite BalrOsQ- We 
found that, in the GGA+U (U I Jff= 0-5 eV) calculations, 
FM alignment of Ir spins was energetically lower than the 
different possible AFM configurations. However experi¬ 
ments indicate a clear AFM transition at 17 K. Following 


the observation of significant orbital moments at both 
Ir sites, we investigated the possible role of SOC in the 
magnetic ground states. For this purpose, we conducted 
GGA+1/+SOC calculations for FM and different AFM 
spin configurations. We found that, with the inclusion of 
SOC, compared to the other possible FM and AFM con¬ 
figurations, the A-type AFM state became the energeti¬ 
cally lowest [| Ea-afm ~Efm |= 2.8 meV/f.u] magnetic 
state. Our findings perfectly agree with previous elec¬ 
tronic structure studies [33| as they also mentioned SOC 
driven stabilisation of the AFM state. In this A-AFM 
configuration, Ir spins are aligned ferromagnetically in 
the ab plane, and were coupled antiferromagnetically in 
the out-of-plane direction. We note that the experimen¬ 
tal band gap (~ 0.3 eV) agrees qualitatively with the 
GGA+ U+SOC ( Ui r ff = 3 eV) AFM band gap (» 0.35 
eV). In Fig. S3 (in the supplementary material), we show 
the GGA+ U and GGA+ U+ SOC band dispersion for the 
FM and AFM configurations. Because of the increase of 
Fermi energy for the AFM state over the FM state with 
the inclusion of SOC, the experimentally observed AFM 
state becomes energetically lower than the FM state. We 
confirmed the stability of the AFM state over FM config¬ 
urations by cross checking our findings across the range of 
f/gjj values from 1 to 5 eV, and we verified the existence 
of the proposed orbital ordering in the presence of AFM 
coupling as well. An interesting point to note is that, to 
obtain the correct magnetic ground state, the inclusion 
of SOC was mandatory, which reflects the importance 
of SOC in this Ir compound. We note a similar mech¬ 
anism was recently suggested for other iridates 33, l35| . 
Therefore this material may be classified as a Coulomb 
enhanced SOC driven AFM insulator. 


It is common believed that strong SOC destroys or¬ 
bital ordering, however, the present material shows a rare 
exception. To understand the possible role of competi¬ 
tion between SOC and structural distortion in deriving 
orbital ordering, we compare our results with the isoelec- 
tronic Ba variant (Ba 2 CeIr 06 ) of the present compound. 
Since Ba 2+ has a larger radius the for the tolerance fac¬ 
tor for the later (t=0.991) is close to 1. Therefore Ba 
form in an ideal cubic crystal structure [36] and the elec¬ 
tronic structure does not allow ai^-e^ splitting required 
for orbital ordering. Because of the smaller radius of 
Sr 2+ , the tolerance factor (t=0.935) of Sr 2 CeIr06 signif¬ 
icantly deviates from 1 and thus a distorted monoclinic 
crystal structure results that allows the orbital order¬ 
ing. From the total energy calculations we found that 
for Ba 2 CeIr 06 , without inclusion of SOC, the monoclinic 
structure is lower (by 138 meV/f.u) than the cubic and 
we found similar orbital ordering as S^CelrOe- However, 
inclusion of SOC cubic structure became « 17 meV/f.u 
lower than the monoclinic. In contrast for the S^CelrOe 
both without and with SOC, the monoclinic structure is 
energetically ~ 747 meV/f.u lower than than the cubic 
one. These analyses gave indirect evidence that, since 
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SOC energy scale remains unchanged between two com¬ 
pounds, in Sr 2 CeIr 06 the structural distortion is strong 
enough to overcome the negative SOC effect and we do 
obtain the orbital ordering state, which appears not to 
be possible for Ba 2 CeIr 06 - 

In conclusion, we conducted both experiments and 
DFT based first-principles calculations to investigate the 
nature of the orbital ordering in the ordered double per- 
ovskite S^CelrOe, with Ir as a single electronically ac¬ 
tive site in the geometrically frustrated fee lattice. We 
found t 2 g band derived anti-ferro type of orbital order¬ 
ing at the Ir site consisting of d yz and d xz orbitals that 
were associated by to alternating Ir sites along the body 
diagonal to the unit cell. Further experiments, such as 
X-ray resonant spectroscopy, would be helpful to directly 
probe our proposed orbital ordering pattern. We found 
that small electron-electron onsite correlation U is also 
relevant in this 5d-Ir system for obtaining the correct 
insulating electronic structure. SOC played a very cru¬ 
cial role in this orbitally active material, for stabilizing 
the experimentally observed AFM ground state in addi¬ 
tion to the small correlation effect, and with a large con¬ 
tribution of orbital magnetic moment greater than the 
spin counterpart. This Ir based double perovskite ma¬ 
terial has a very striking feature in which both electron- 
electron interaction and SOC energy scales are important 
for understanding the electronic and magnetic behaviour. 
More detailed experiments and theoretical analysis are 
required to explore this critical energy balance in this 
interesting material. 
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